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FOREWORD 


This  interim  technical  report  was  prepared  by  the 
Aerophysics  Laboratory  of  the  Department  of  Aeronautics  and 
Astronautics,  Massachusetts  Institute  of  Technology,  Cam¬ 
bridge,  Massachusetts  on  Contract  AF  33(657)-7975  for  the 
Aeronautical  Research  Laboratories,  Office  of  Aerospace 
Research,  United  States  Air  Force.  The  research  reported 

m 

herein  was  accomplished  on  Task  7065-01,  "Fluid  Dynamics 
Facilities  Research"  of  Project  7065,  "Aerospace  Simula¬ 
tion  Techniques  Research"  under  the  technical  cognizance  of 
Mr.  Robert  G.  Dunn  of  the  Fluid  Dynamics  Facilities  Labo¬ 
ratory  of  ARL. 


ABSTRACT 


One  possible  technique  for  accelerating  a  plasma  to 
a  high  Mach  number  is  that  of  passing  a  magnetic  field 
through  the  plasma.  This  technique,  called  the  traveling 
wave  pump,  is  characterized  by  a  complicated  set  of  differ¬ 
ential  equations.  These  equations  have  been  approximated 
by  a  one-dimensional  steady  state  form.  Several  approxi¬ 
mate  integrals  are  found  for  the  one-dimensional  equations. 
The  results  indicate  that  the  viscosity  losses  in  the  system 
are  not  excessive.  The  results  also  indicate  that  the  inlet 
velocity  profile  should  be  as  uniform  as  possible  and  that 
the  magnitude  of  the  inlet  velocity  should  exceed  /-//>  times 
the  wave  speed  by  an  appreciable  amount. 
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CHAPTER  I 


ANALYSIS 


The  complexity  of  the  flow  in  the  TWP  accelerator  seems  to 
require  a  simpler  procedure  for  a  working  approximation.  Since  the 
TWP  is  an  A.C.  apparatus  the  first  step  is  that  of  reducing  the  equations 
to  the  form  of  steady  flow.  The  one-fluid  model  will  be  used.  The  next 
step  is  that  of  reducing  the  equations  to  one-dimensional  flow.  These 
steps  which  give  simpler,  though  more  approximate  set  of  equations 
are  outlined  below. 

Consider  first  the  continuity  of  mass, 

&  -  </r  r  £0  9)  -  O  0) 


where 

/O  -  mass  density 
p  =  mass  velocity 
t  =  time. 

e 

If  this  equation  is  integrated  over  a  volume  (this  volume  is  not  time 
dependent)  and  the  divergence  is  converted  to  a  surface  integral,  one 
obtains 

f  4^—  cf\/  t  f *  ds  *  O  (Z) 

ypf  s 
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Equation  (2)  !■  then  integrated  over  a  period,  i.e., 

/  e/t  fj  "df"  c/v)  +  / dt  f av  *  c/s  *o  (S) 

t  »**/  St  S  ' 


and  if 

/  0  "I '/  fr  ", / 

*»/  V  Mb/  't+T 

then  thia  term  vanishes.  The  second  term  is  evaluated  over  a  stream 
tube  of  cross  section  P  fz)  lying  between  Z  and  Z  +  dz 
Thus  the  continuity  reduces  to 


yOiz  S?  *  Cons/-  or 

(■**.) 

* 

* 

* 

ii 

* 

(*b) 

This  volume  corresponds  to  the  control  volume  used  by  Shapiro. 

Similar  treatment  of  the  momentum  equation  results  in  the 
average  one  dimensional  equation, 


dP 

P  P4  2 


O 


where 


M  =  local  Mach  number 


P  =  pressure 
^  =  ratio  of  specific  heats 


2 


•tl  -  total  body  force  acting  on  the  control  volume  in 
direction  opposed  to  U. 


o 

o 

f 


=  hydraulic  diameter 

.  area  of  duct 

=  4  - 

perimeter 

=  friction  factor 

=  rw/t  i  *  1  T'ou‘ 


The  energy  equation  is  also  averaged  over  the  period  and  gives 


7-/ 

2 


where 


4* 

dQ 


-  specific  heat  at  constant  pressure 
=  work  and  heat  added  to  the  contents  of  the 
control  volume  per  unit  mass  of  contents 


The  state  equation  is 


-  _f/0  .  d/>  dT 
M  ~  fi  +  4  *  T 


The  next  step  is  that  of  evaluating  the  averaging  integrals,  which 
can  be  calculated  from  the  results  given  in  Ref.  2  after  computing  the 
electric  and  magnetic  fiefds.  When  the  magnetic  Reynolds  number  per 
unit  length  based  on  the  slip  velocity  (  CC~U.)  ),  is  small, 

the  attenuation  of  the  coil  is  small  and  (c-u)  is  small.  The 

arguments  of  the  Bessel  functions  describing  the  fields  below  are  real, 


3 


with  M  *  ^ 


.  The  fields  are  given  below  for  this  case: 


sr  It  fxr)  Gfx*c)  S/nll 

e2  =  yuASZ  ZQ  f*r)  $fX8j  Cos  a 
=  -CyjLAtl  f(**c)  If  OCr)  sin  (X 
£  =  -01/aA/Z  l/*r)C(**0)  (e-vz) 


The  force  per  unit  volume  is 

F  *  3  *  3rJ0  JQg  6rTs 

Fr^-<r  f^u  XI)*  0  *(**c)  Ia  (*r)  If  filr)  sin  A  Cos  ZX  (C- 1£) 
Fz  «  <T(yu/VI)*  Z,  *(Kr)G*f**J  sin  *A  (C-  Vx) 


The  total  electrical  energy  added  per  unit  volume,  per  unit  time  (i.e.,  the 
total  power  addition)  is 


PT  -  t-J  =  <r CmA/T)J  Tt  *f/(r)6*0(/><>)  sio*Si  C/C-Vx) 


where 

K  *  i/c  ,  A  •  K  fz -  ct ) 


4 


and 


)  * 


_ K  f**t) _ 

^  (**.)  it  f**t )  -Ka  (**,)  zt  (kk,) 


where  I and  JCf^^  are  Bessel  functions  of  the  first  and  second 
kind,  respectively,  of  order  A/  having  imaginary  argument. 


For  KR. 


2+Oc*)2  ln(KR) 


X  / 


The  force  on  the  control  volume  of  length  Z  then  becomes 

z+A  t+r  4 

£  =  dxT f*nJG  *r*V  'C**1*  *h,  *•”**(*- ct)  rdr 

*■  t  • 

*  -<rfiuA/z)*~  •  <Kr)(C-vJr)rd/* 

s  T  <tC+l#Z)2  ~I0  f*Q ]  times 

G  *(*%)  m  £c-txu]  z 

Zn  f  \  *(Kr)  Vt  (r)rdr 

n*e*rx}(KV  -I,  f*Vzx  ,«v] 


Of  is  a  weighting  function  which  adjusts  (/•  ,  the  mass  average 

velocity  to  account  for  the  variation  in  magnetic  field  across  the  tube. 
In  the  case  of  slug  velocity  =1. 


5 


Let 


Ft*:*)  -  [l*  t kb)  -I0  r**)  i  fKR)  ] 


Thin  as  tt*  — ►  oo  J  Ft  A'*) 


Si • 


J## 


frrt*#)* 


and 


cfZ  =  -*^*x)*  Q *tk'/y  rr#fFtx*J(c-ocu.)  dz 

and 

^  ’  T  I/z  [rf^  (f  H  *)]  J* 

=  -f~-  Q  (K*t)  r <■*#.)  Cfc-tu) 


A  one  dimensional  form  has  also  been  deduced  by  Williams  in  Ref.  3. 
However  he  fails  to  account  for 

Equations  (4),  (5),  (6),  and  (7)  describe  the  process  to  be 
analyzed.  Since  these  equations  constitute  a  non  linear  system,  it 
seems  most  practical  to  follow  Dahlberg's  suggestion  in  Ref.  4  and  study 
the  phase  plane  in  search  of  closed  form  solutions. 

The  momentum  equation  then  becomes: 

'£'T*r  +  i'M*p{2-L+1Lg),0  fg) 


6 


Using  ^ 


* 


7  Pm* 


/’“Jr *57  .  ir^1-  F(K*j(}'ot*.uc-<>'u)-Zfiu.*-% 


or 


^  V*^;  fC-oeu)  -'*p**'£‘  '‘/0UJ~  (& 


The  energy  equation  (6),  neglecting  the  ionization  term, 

c/  UL  ^  Zd O* 

equation  of  state  and  expanding  -  ——  =  ■■■^  1 


using  the 
,  becomes 


dT  p  C+ 

r  /Oje  P 


u.du. 


From  (4b)  and  (7) 


dT  dP  do  dP  ^  du.  .  dp 

~  -  —  -ir*  —  *—*— 


so 


dz  ' 


p  du  Pu.  dp 
*  rd*  f~dT 


*• 


7-/ 

P 


/>u 


2  du 
dz 


] 
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leading  to 


du  t/P 

dz  y  /f  dz 


c*  Pu  c/0 
X  *  dz 


j  *Cu»r)* 


r(x*A  )G*f*Xs)  CfC-  aca) 


Solving  for 


du. 

dz. 


dz 


_  <r(#/*T)z 
+ 


$*(/<*.  >  c(c-  *u)  ^ 


/>!**• 


P 


Substituting  in  ('))  and  collecting  terms 


dP 
dz  " 


-  [/OU.**^-  p]  F(K*9)G*(*eJ(C-<Xu) 

+/OUCT  FfxR* )  <f  *(*/?•  )  c (C -<XU)J 

[■&  -/]  />“*  ~(^k)p 


0?) 


8 


Examination  of  these  two  equations  indicates  that 

a.  there  is  a  critical  line  corresponding  to 


or  since  p  ,  and  P  are  not  generally  zero  this 
line  corresponds  to  Aff  ■  / 
b.  there  are  singular  points  where 
occurs  if 


*7 


This 


1. 


U  m 


and 


ii. 


*  “ 


cr*  _  -  . 
Zz  “  Z/°' 


O 


When  this  singular  point  is  located  within  the  region  of  interest  it  offers 
a  possible  termination  for  the  process. 


9 


CHAPTER  II 


DISCU  SSI  ON 


First  consider  the  equations  (II)  and  (12)  in  the  physical  plane. 
These  two  equations  can  be  integrated  numerically  from  given  conditions 
at  the  inlet,  to  the  outlet  of  a  TWP  device  once  its  geometry  has  been 
specified. 

In  order  to  find  the  circumstances  which  will  lead  to  an  increase 
in  Mach  number,  consider  the  following  expression  from  Ref.  1  which 
follows  from  the  momentum  state  continuity  and  energy  equations. 


20  + 


r-t 


M*) 


*  /-A4* 


c/Q 


rht*Q+^irf*z) 

f-M* 


/ 

2 


This  can  be  rewritten  as: 


X  /*  -IT 

to  ctx.  ~  /-M* 


i**££±  dt>  l  *  l iSLal  / 

/-M*  a'*  iCfiT  /-M*  e/x 
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Consider  the  case  when  the  friction  cancels  the  area  change. 
This  case  exhibits  the  effects  of  the  magnetic  power  addition  only. 

dA4 


For 


dz 


to  be  positive 


~  £. 

2  2  Ota 


03a.) 


dM 

For  this  case  ■»  a  O  lines  for  constant  values  of  IT  are  shown 
in  Fig.  1.  These  lines  separate  the  plane  into  <  Q  and 

regions.  Also  for  this  case  if  —  is  also  to  be  positive,  U.IC 

must  be  greater  than  -*■ — — 


This  constraint  is  considerably  less  stringent  than  the  previous 

du.  d  M 

one;(  can  be  positive  when  51  is  0  or  negative)  hence  it 

will  not  be  plotted.  dz 

Oahlberg  (Ref.  4)  has  shown  that  Eqs.  (11)  and  ( 1  2)  can  be  put 
into  a  symmetrical  form  by  the  proper  choice  of  nondimensional 
variables.  ^  « 

Choose  VL  a  — -  and  ^  *  jQ  U.C  as  new  variable®  and  f/fj 

express  the  equations  in  terms  of  two  nondimensional  parameters: 


>  *  5*75 


and  ^  » 


#e<r 


JL 

o 


0+a) 


and  the  characteristic  length  — 


/OS4. 


For  the  case  of  constant  phase  velocity  £*  and  gas  properties  the 
equations  become 


du*  ) +u*fdu*-y,p+) 

A.  "  s  —  J  ^ 

/oi 


(>s) 


11 


e/P 

dz 


(oeu  *-/)(P  *-/)  +{P*+u*)(4**~  "WP  V 

p+-r-/  +  -!£ )  u* 


For  the  case  with  cylindrical  symmetry  is  the  equivalent 

average  field. 


above, 


Tn  terms  of  these  variables  at  the  singular  point  mentioned 
^  a  /I  ,  i.e.,  the  nondimension&l  pressure 

dz  dz 

and  velocity  become  constant  giving  in  a  sense  a  "fully  developed,” 

stationary  solution,  although,  of  course,  the  dimensional  pressure  must 

decrease  along  the  duct  to  balance  the  friction  under  this  condition.  The 

nondimensional  variables  thus  give  a  much  better  description  of  the 

termination  point  of  the  process. 

This  .singular  point  toward  which  the  magnetic  pumping  drives 

the  flow,  must  be  located  in  the  supersonic  part  of  the  phase  plane  and 

the  friction  must  not  be  too  high  if  supersonic  acceleration  is  to  occur. 

•  - 

These  limits  are  expressed  by  regions  t  and  it  of  the  ff  •  ^ 
plot  of  n.ef.  4  and  impose  the  constraints: 


a. 


_ /_  e/P 

2  *P  c/z 


6. 


***** 

*/>u  rfr-t) 


* 


JP 

O 


('*) 
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CHAPTER  III 


SOLUTION  IN  THE  PHASE  PLANE 

Equations  (15)  and  (16)  are  identical  to  those  of  Dahlberg  (Ref.  4) 
when  the  variation  of  the  field  with  radius  is  accounted  for  through 
and  OC  the  velocity  averaging  factor.  We  follow  his  procedure 
dividing  (16)  by  (15)  to  obtain 


d? 

du 


('ecu-  f)(P  +  U-0 
(ccu-O(u-JL)  /  uftu-yrp) 


The  *’  s  have  been  omitted;  however  the  nondimensional  variables 
(defined  in  (13)  and  (14))  arc  still  being  used. 

Following  Lighthill  (Ref.  5)  let 

P'Potr)  +  .  .  . 

and  (/Q) 

*  *  X  +  Vtty  (*)  +  ..  . 
efP 

dP  27 

du.  du 

dx 


c/u 

dx 


=  /  v*  •  • 


dP 

dx 


a 


r.'+vr.'- 


•  • 


09) 
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(17)  can  then  be  re-written 

)  +u  (0*  “  V Pjj  =  M  jfru-/)fP+u-t)  + 

+  SP+u) (Pu  -  VP)]  (26) 


Substituting  (18)  and  (19)  into  (20)  and  equating  powers  of  V , 
the  nondimensional  area  change  parameter,  we  obtain  the  O  order 
equation  for  P  *  * 


_p,n  «*X-/)(P0  +*-')  +  P*(P0  +*) 

°  (ou -/)(*- -£  )  +  px* 

>» 

and  the  two  first  order  equations: 


( 2/ ) 


p'_.  p  (<**-«)  X-! _ 

'  '  «XX-/)  (X- £-)  +  &  ** 

cf> 

and 

u'  +-fU)Ut  * f(x) 

where 

o<rPm  +*-/)+ foe*  *  r- /+?p0 

(X)  "  fax-/)  (P0  t  x-t)  +  j  z  (/»  f-  z) 

(<x+p)x  +fa-/)x -(/+  %£) 
fax-/J(x-£)  +  pr* 


and 


7>  (P0  +  x) 


9CX)  *  fa*-/)  rp0  +x-/j  +*x(p0  +x) 

_ _ 

(* x-/)(y-£  )+px* 

CP 


(22) 


(23) 


~\ 
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Equation  (21)  is  a  first  order,  linear  differential  equation  for 
which  can  be  reduced  to  quadratures  and  integrated  in  the  case  when  ^ 
is  small.  /^cO  is  a  regular  solution  of  (22)  which  is  chosen  to  sat- 
isfy  the  initial  conditions  that  at  it  t  PK  ^  ^ Pf  s  P-'jf  f 

P,*o. 

After  solution  of  (21),  (  23)  can  be  reduced  to  quadratures,  but  it 
involves  a  great  many  terms  and  an  approximation  to  the  logarithmic 
function  even  for  small  and 

a.  The  Solution  for 

Collecting  the  terms  in  (  23)  we  obtain 


„  («**)*•/ _  ..  fattix*.  ifftt*/)*/  .  . 


The  solution  of  this  equation  is 


M  •  exp  C-TMjfe,  Tfrj  («!*)**- *(*<■»  *  > 

V  J  (,*%)+£  J 


(ZS) 


where 


I(t) .  / - _ 


c/jc 
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This  can  be  integrated*  to  give 


•TM  *•{/»  p*+f)**-*0+jr*  >+tr  J 


using  this  result  (25)  becomes 


(Zt) 


where  the  following  abbreviations  have  been  used. 


*  For  the  Logarithmic  function  to  be  the  solution 
is  required.  Rewriting  this  condition  as 


it  can  be  seen  that  with 

^  *  ^pr  /±r*/.*7 

thi3  condition  is  nearly  always  met.  For  example,  the  limit  is  •  ZZ6 . 
For  the  worst  condition  of  flf  7**  /• 
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A*  -ft*  £r) 


°-e?> 


4-B  «<*/"/ 'cc+4)>t 


9> 


c+ 


The  second  term  of  (26)  cannot  be  integrated  in  general;  but  when 
small  (see  Appendix  II),  the  expression  can  be  integrated  directly  since 

/??—►/,  n  -+-0. 

F  or  and  absorbing  a  numerical  factor  in  the  constant  Cf  . 


P. 


(a.  +f)r-  24  -  f.(& 

ei*  ,.z* 


is 
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ot+4 


(OL++)X- 


/ 


(eL  +  (*+t)x)*'  - 


«?**«+/ 

(«**)* 


(a.+  (b*4)x)  * 


&*+ou(ot+i>+ot+4 
(ot+4) 1 


•In  (-A-fm-4) g)x 


(27) 


where  fl  *  " 

(*+<*+*)*) 


/  ^  ■■  ^  . 

appears  to  be  negative. 


For  any  physical  situation 


b.  Higher  Order  Terms 

Since  the  substitution  of  (27)  in  (23)  leads  to  a  quadrature  which 
offers  little  advantage  over  direct  numerical  integration  in  the  physical 
plane  for  a  specific  geometry  device,  further  attention  is  directed  only 
to  the  simple  case  of  constant  area  which  is  integrable  directly, 
a)  Constant  area  <JP  a  O)  no  friction  (+•<>) 

When  )►»(),  tfi  /  #  P  •  P9  the  general  phase 
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plane  trajectories  are  shown  in  Fig.  4a  of  Ref.  4.  With  these  simplifica¬ 
tions 
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where 


C ^  is  of  course  determined  by  the  condition  that  at  the  inlet  1*0  , 

P  *  PfnUt  ,  *  m  “intmt 

The  rest  of  the  right  hand  sidej  evaluated  at  U  •  &tn  • 

b)  Y  *0  (constant  area)  ^  «<  /  (small  friction)  the 
general  phase  plane  for  this  case  is  shown  in  Fig.  4b  of  Ref.  4.  Super¬ 
sonic  acceleration  in  this  case  is  confined  to  every  small  region  of  the 
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phaac  plane  and  can  occur  only  when  <t>  <  !/¥?(*-!)- 
Then 

/ 

- 


(a) 


(a  +  4>)u  - 


OR 


[c' 


+  —  ( a  +  (&+4)u) 


a  +  t 

a1  ¥  a  (<t+o  +  (a+t) 


a  +  d> 


Ln  -[a  +  (a+<f>)u\  J 


where 


C!  =  ^//l 


* 


tf/? 


-cl  c/> 


+  2(d+<i»  (<i  *  uih)Z  ~  (a  +(<*+&  um) 


a+<t> 

For  this  case  (15)  gives 

dz  - 

3 
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du 


(30) 


(31) 


neglecting  powers  of  <fc  greater  than  the  first  (30)  becomes: 


_ / _  r  (au  -D* +2(au-l)4>6 

(a  +  4>)  u  -  —  -  j>6  L*  2(a+<j>) 

cp 

a-i+zts 

+  - - -  (ctu-/ 1  4>6) 

a+  <f> 


~  ~atf  (<*6  +  6  +  0  Ln  (/-CLU) 1  (32) 

where 


J2- 


integrating  (3  1 )  with  (32)  we  obtain 


z 


It  *  c'h  ~  7<a7ft 


a-/  +  2+6  ,  t«t6  +  6+/)  r' 

+  a+t -  - VT*  J'J  (33) 


where  all  the  integrals  can  be  found  explicitly  exeept  Xf  which  requires 
the  approximation  L  n  ( f  -  Ct  U)  m  -  C&  U  +  €  (CL  u)*)  where  € 

is  taken  to  give  the  best  approximation  to  the  L  n  over  the  range  of  Ot<x 
of  interest  (sec  Appendix  II).  For  Ctm  ^  . S  €  *  //<2  gives 
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about  4  per  cent  error  under  the  worst  condition  CdU*.S)  •  Since  this 

term  is  of  order  <f>  the  error  introduced  in  Z  should  be  small.  The 
above  integrals  are  as  follows: 


/•  d  u 

r 

(a +  4>)u*  - 

[<fa  *  <P)  u  -  ( /  - 


du 


«t +4) 


AiS 

_  1  . 

A2S 

i 

8* 

[  ] 

“T 

[  ] 

V 

[  ] 

du 


(a+$) 


Az*  .  *4 


t  ]  [  ]  t  )\ 

*£ 


du 


J  It  ]  t  ] 


[  i 


du 
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where 


q[|Sf  0-2*)+  <lt<p^U-la)+  Z-^-(ai-l)fZa-aS)^ 

[fjf  ~  l  +  2*6](a++)t 
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a  (/-a)  y-  <f>  \aa26  -  za(S\ 

[f^f-  ~  1  +  2  <fi)z 
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A2s 

d  -  a2 ( / -&  +*)  +p(a -za*€6) 
(§£•  -1+246)*  (a+?)2 

with  the  above  substitutions  (33)  becomes: 
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{3}  Ln^(ct+4)u  - (I  -  <f>6) ] 


^  0  (ol+4>)3 


(3+) 


where  [  everything  else  on  the  right  hand  sidej  evaluated 

at  U  *  U  in/ 9 1  The  {  ) 's  are  constants  involving  Jr  t  CL  ,  $ 

and  €  the  Logarithmic  approximation,  specifically 


{/}  =  -)  +  t  at-»a*0-aS+6) 


-Y<*6-a*S  - 1  a* (j-j* (a6  +  6  +0 


d*R  1 

+  a<5(*(t  +  at*  +!)-<*  +2  ~c~~  (a-f-t)j 


(3}  =  StiLi*L  +  ^  |^ar<f  (a1  +n-dl6(aS  +  8+D 

+  6 a*  (c t6  +  6  +  /)j 
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CHAPTER  IV 


CONCLUSION 

The  expressions  found  by  the  method  of  Lighthill  for  P(u)  and  it  (z) 
for  one  dimensional  compressible  flow  at  constant  area  should  facili¬ 
tate  the  interpretation  of  the  experiments  with  the  constant  area  TWP 
with  small  friction.  Although  this  method  produced  a  quadrature  for 
the  first  order  component  of  velocity  due  to  area  change,  it  seems 
doubtful  that  evaluating  this  expression  would  result  in  any  simplifi¬ 
cation  over  direct  numerical  solution  for  a  specific  case. 

The  results  also  indicate  the  regions  of  operation  in  the 
local  Mach  number  plane.  This  plane  illustrates  the  conditions  that 
must  be  satisfied  to  increase  both  Mach  number  and  velocity  in  this 
kind  of  process.  This  curve  seems  to  offer  the  suggestion  that, 
according  to  the  one- dimensional  flow  it  will  be  necessary  to 
increase  the  phase  velocity  of  the  wave  as  it  travels  down  the  accel¬ 
erator.  The  one-dimensional  approximation  also  indicates  that  the 
performance  is  improved  by  making  the  inlet  conditions  of  the  gas 
as  uniform  as  possible. 

The  calculations  in  the  appendix  indicate  that  in  terms  of  Dahl- 
berg's  nondimens ional  friction,  the  friction  will  be  small  in  devices 
operating  at  low  pressures.  The  advantage  of  maintaining  a  disequi¬ 
librium  be.tween  electron  temperature  and  gas  temperature  is  apparent, 
since  elevated  electron  temperature  gives  increased  conductivity  with¬ 
out  the  increase  in  viscosity  which  accompanies  increased  gas  tempera¬ 
ture.  Some  data  is  also  presented  in  Appendix  I  which  serves  to  compare 
various  gases  in  terms  of  the  friction  parameter  ^ 
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Figure  1.  Conditions  for  local  change  in  Mach  number 


APPENDIX  1 


THE  NONDIMEN  SIGNAL  FRICTION  FACTOR 


^  -f 

Dahlberg  (Ref.  4)  defines  the  friction  parameter  P  -  J  ■  ■ 

in  terms  of  the  friction  factor  T  .  In  terms  of  the  skin  friction  coeffi- 


»  j  _  /»•»  v  ^ 

cient,  C^.  ,  this  becomes  0  can  be  approximated 

as  the  laminar  or  turbulent  skin  friction  coefficient  depending  on  the 


magnitude  of  the  Reynolds  number,  s 
gas  law  it  follows  that 


-  /*««  O 


From  the  perfect 


x,  /jSlT)  =  ~m.  J&~>  *.  0,66 

r  \  AT*  J  roc) 


where  P  is  in  mm  of  mercury  T  is  in  °K  and  p  is  in  AfAf/ AT* 
is  the  molecular  weight.  Using  the  simplified  kinetic  theory  of  Hirsch- 
felder,  Curtis  and  Bird.  (Ref.  6) 


-  .. 

2.67*  >0  ~S7~ 


where  is  the  equiv¬ 


alent  rigid  sphere  molecular  diameter  in  Angstrom  units.  Substituting 
into  the  Reynolds  number  wc  obtain 


uO  — *  6.02  *,o * 


for  T  at  5000  *K  and  1  mm  pressure  with 

D*./m  a  «  / o+m/sec  =  5/ 

For  Reynolds  numbers  of  this  order  the  laminar  skin  friction  will  certainly 
be  applicable  in  which  case 


With  the  above  relations  the  friction  parameter  becomes 
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j.  /.  J*r/o  -*  Afr  Plmm>  * 

~  t'*-*,*,  V5^  v<V'" 

for  nitrogen  at  5000  ®K  and  lmm  pressure  with  /  weber/meter 
and  <T  =  lOOmho/fTt  j  »  *06  X  /0  ~S 

In  Table  I  the  appropriate  constants  for  some  commonly  used 
gases  are  given  to  facilitate  a  comparison  of  their  performance  in 
the  TWP. 


TABLE  I 


Proportional  to 


UW 

•V*L 

m 

fas 

^4° 

•  V 

tFZ** 

Cross  Section 
(Angstroms ) 

Molecular 

Weight 

(*) 

(Re) 

N* 

3.75 

28 

.376 

74.  5 

3.  25 

He 

2.  18 

4 

.42 

9-  5 

1.3 

A 

3.64 

40 

.48 

84 

4.4 

h2 

2.92 

2 

.  166 

12.  1 

.57 

Air 

3.7 

29 

.394 

74 

3.4 

Hg 

3.  56 

200 

1.  14 

180 

15.  1 

N 

1.68 

14 

1.33 

10.5 

4.3 

co2 

3.89 

44 

.44 

100 

4.37 

ch4 

3.79 

16 

.278 

57.5 

2.  1 

From  the  above  table  it  appears  that  hydrogen  has  by  far  the  best 
characteristics  as  far  as  the  frictional  losses  are  concerned.  For  most 
of  the  cases  of  interest  the  small  ^  assumption  appears  to  be  justified. 


**  This  equation  can  also  be  written  in  terms  of  the  Mach  number, 
u/a  .  It  becomes 

r  <r£M^fo  <^s 
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APPENDIX  II 


APPROXIMATION  TO  THE  LOGARITHMIC  FUNCTION 

In  order  to  approximate  ,  consider  the  first  two  terms  in 
the  series  for  *  -  [/Qru )  *  ^  ft* U)  X.  {QfU)S.  .  .] 

In  order  that  the  integral  be  expandable  in  partial  fractions  the  numer¬ 
ator  must  be  of  order  2  or  less.  We  choose  the  approximation 
£/>  (t  -  ctu)  V  -  [ocu+f  1»u)*J  where  €  is  matched  at  a 

point  near  the  average  o {OtU.  over  the  acceleration  range  of  interest. 
For  OtU  <  S €  *  gives  a  good  approximation,  the  range 

46  - 

of  validity  of  a  particular  value  getting  worse  as  OtU  s  /  is  ap¬ 
proached.  Values  of  £  for  a  perfect  fit  at  a  given  value  of  OtU.  are 
given  below. 

v  2 

OtU  *  for  lh(i-ccu)*-<xu.-€(oeu) 

.4  .69 

.  5  .77 

.6  .89 

.7  1 . 02 

.8  1.25 

.9  1.48 

.95  2.28 
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APPENDIX  III 


THE  FIELD-VELOCITY  WEIGHTING  FUNCTION 


Because  the  fields  in  the  TWP  vary  so  strongly  with  radius,  if 
is  large  the  field-velocity  weighting  function,  0£  ,  must  be  included.  In 
the  text  OC  ia  defined  by 

~f-  -  ■C  yfX'dr 

h, *(*».) -T,(**C)ZX  («*.>] 

To  predict  the  performance  of  a  TWP  device  by  means  of  the  one-dimen¬ 
sional  approximate  analysis,  a  knowledge  of  the  behavior  of  Oe  for 
different  velocity  profiles,  and  various  magnitudes  of  MR0  is  required. 

Assuming  the  density  is  constant  across  the  cross  section,  46  , 
the  mass  average  velocity,  is  V (t) 

O  % 

As  a  simple  representation  of  a  general  velocity  profile  assume 

V(r)  =  u6  o±r  ±re 


as  represented  in  Fig.  A-l.  Then 


and 


a 


-z  'W  ] 
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ifg  x,  r*jr  +*Jr  * 


The  last  term  in  the  {}  requires  a  numerical  integration  or 
an  approximation  to  I  For  small  argument  (**„  “O 

Te(**.)z /  +  (■*£)* 

I.  t**4)  3  ~ 


jr  (**J  Z  -£  (^-)* 


and  OC 


TJi  {*(>-#)')} 


Some  values  of  Oc  | ft. 


for  various  assumed  values  of 
r. 


ffo  are  plotted  in  Fig.  A-l.  Note  that  as  -  ►/  oc  must  always 
approach  /  from  below  regardless  of  the  value  of  ,  and  at  5 

this  approach  is  very  rapid  from  O  in  the  immediate  vicinity  of  -A  ; / 
Because  of  the  form  of  the  expression  for  OC  the  asymptotic  approxi¬ 
mation  to  JTt  does  not  give  the  right  limit  as  ^ ■»/. 

Figure  A-l  indicates  that  for  most  applications  there  would  be 
little  advantage  in  operating  above  ^  f  unless  a  very  nearly 

slug  velocity  profile  was  present. 
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Figure  A-l.  Illustrations  of  effect  of  nonuniform  profiles 


